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ABSTRACT 

We investigate the effects of H-dibaryons on the structure of neutron stars. We find 
that H-particles could be present in neutron stars for a wide range of dibaryon masses. The 
appearance of dibaryons softens the equation of state, lowers the maximum neutron star 
mass, and affects the transport properties of dense matter. We constrain the parameter 
space for dibaryons by requiring that a 1.44 M 0 neutron star be gravitationally stable. 
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The existence of long-lived (or even stable) six-quark states with baryon number 2 was 
first proposed by Jaffe 1 on the basis of bag model considerations; experimental searches 
for them are being actively pursued,* with one recent report of a positive detection. 3 Of 
the possible dibaryon states, the H-particle, with the quark content of two A’s (uuddss) 
in a flavor-singlet spin-zero state ( J p = 0+), is the best candidate for metastability. Bag, 
Skyrme, quark cluster, and lattice QCD model calculations predict a mass for the H- 
dibaryon in the range m B ~ 2.1 - 2.24 GeV, with most models indicating a value 20 - 80 

MeV below that of two A’s (2m A = 2.231 GeV). 4 ^ ttis case ’ the H'P article would be 
stable with respect to the strong force, decaying only via weak interactions. 

Dibaryons will be present in neutron stars whenever the chemical potential of neutron 
star matter is greater than the effective mass of the dibaryon. Since they are bosons, like 
pions and kaons, dibaryons can, in principle, condense in the cores of neutron stars, but 
with one important difference: dibaryons carry baryon number. If dibaryon interactions 
are neglected, a Bose condensate of dibaryons would cause the collapse of neutron stars 
with central densities higher than the minimum density for the appearance of dibaryons. 5 
At some level, however, such interactions must be present, and such an instability will 
not be realized. At very short distances, the interactions among the underlying quark 
components of the dibaryons will give rise to an effective repulsive interaction between 
dibaryons. Even a very small self-coupling can stabilize the system against collapse. If 
the interactions are relatively weak, however, the equation of state of dense matter will 
be softened by the presence of dibaryons. 

At a critical baryon number density, n B > »J, dibaryons begin to appear inside neu- 
tron stars. We study this transition by considering a mixture of neutrons and dibaryons in 
chemical equilibrium. (For self-consistency, we should also include protons, As, and other 
light hyperons, but this does not qualitatively change the results for the parameter range 
of interest. 8 ) In equilibrium, the chemical potentials of dibaryons (fi H ) and neutrons (/i n ) 
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are related by 


Ha = 2/*„. W 

Of the many candidate equations of state for neutron matter, we choose the medium- 
stiffness Bethe- Johnson 7 (BJl) model as a representative and analytically tractable ex- 
ample. (Clearly, it would be of interest to make a systematic study with different nuclear 
equations of state; here, we have chosen a familiar model to qualitatively illustrate the 
effects of dibaryons.) In its polytropic incarnation, we can write the energy per neutron 
as 

e n = — = m n + 236n“ MeV, (2) 

n n 

where the neutron number density n„ is given in fm -3 , m n is the neutron mass in MeV, 
and a = 1.54. The neutron pressure and chemical potential can be written as: 

Vn = nip- = 363.44n“ +1 MeV fm" 3 , (3) 

dn n 

Hn = m n + 599.44n“ MeV. (4) 

The adiabatic index is then T n = a + 1 — 2.54. 

We shall assume that the H-dibaryon is spatially compact, with a size comparable to 
that of nucleons. This assumption is borne out by bag model estimates 4 and corroborated 
by a recent study® of the quark structure of the H-partide. Then, at the nucleon level, 
in nuclear mean field theory, the H-particle can be modelled as a complex scalar field, (f>, 
with effective Lagrangian 9 : 

L t = - l m , w 

> 3 describes higher order multi-body interactions, and Lint includes the coupling 
to other hadrons. In mean field theory, the primary effect of the H interaction with 
nucleons is to induce a shift m H -* m* B in the effective H-mass, due to the mean field 
of the neutrons. Although density dependent, to first approximation we can absorb this 
effect in a mean renormalized mass for H; for simplicity, we will denote this effective 
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mass by m u hereafter. As a first approach to modelling the dibaryons that is simple but 
qualitatively reliable, we will neglect the higher order self-interactions between dibaryons. 
(These terms may be important at high densities, but in that regime the description of 

the EOS in terms of nucleons is suspect at any rate.) 

We can estimate the dibaryon self-coupling A by using a variant of the P-matrix 
formalism of Jaffe and Low, 10 applied to diquarks by Donoghue and Sateesh. 11 If we 
consider two dibaryons brought together, the energy of the resulting configuration can be 
written as E 2 = 2m H + A% A), where r is their separation. As r -» Ra, E, -* M it where 
M 4 and Ra are the mass and radius of the 12-quark baryon number-4 state. The mass 
Ma can be estimated from the bag model, 4 M 4 ^ 4.7 GeV. Calculating the interaction 
energy for the two-dibaryon state confined inside a cavity of radius R, we find that 11 
E int (R, A) = 0.34A/4E 3 H 3 , where E* = p 3 + and the momentum p = v/R from the 
vanishing of <f> at the cavity boundary. Taking the cavity to be the size of the baryon 
number-4 state, R = Ra, and equating AE(Ra ,A) = E int {RA, A), we obtain 

A = 11.8f$(Af 4 - 2 m H ) ^ • ( 6 ) 

For the range Ra = 1.1 - 1.4 fm and m H = 2 - 2.2 GeV, we estimate A ~ (3.4 - 13) x 10 3 . 
Given the crudeness of the approximations involved, this estimate should be considered 
a guide to the expected order of magnitude for the coupling. 

Since the self-coupling A lies comfortably above a critical value A* = 4ir 
~ 4 x lO’ 37 , we can use the equation of state for self-interacting bosons derived by 
Colpi, Shapiro, and Wasserman 13 to excellent approximation (corrections are of order 
C?(A* /A) ~ 10 -40 ). For a spherically symmetric scalar field with time-independent energy- 
momentum tensor, the lowest energy solution is of the form 

^(r,t) = $(r)e- i “‘ , (7) 

where $(r) is a real function of the radial coordinate. In the limit A > A*, the solution 
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( 8 ) 


to the scalar field equation of motion to lowest order in A*/A is just 

/ 2 \ 1/3 

$ = ™ [^L _i) for r<R, 

A 1 / 3 \m 3 y 

where R denotes the radius of the configuration. 13 The dibaryon equation of state is 
obtained by substituting this solution into the scalar field energy-momentum tensor, 

T£ = (dr<fj dvt + d 9 <t>d v ^ - i d x <$ d v <j> + m 2 H (f^<j> + -A(^<£) 2 J (9) 

and the H-particle current 

J» = l g d v <t> - . (10) 

Defining the dimensionless parameter z = u; 3 /m 3 , the dibaryon energy density and pres- 
sure are then given by: 

Pb — Po(-® — l)(3z 4- 1)> (11) 

Ph = Po(z — 1 )*» ( 1 ^) 


where 


4 / . » „„„ MeV f mg 

Po = m B j 4A = 523^3- ^joOOMeV 


')'(?)■ 


The parameter z is related to the baryon number density of dibaryons by: 


ns — 2ng = -^—(z — 1)n/z. 

TTlff 


(13) 


(14) 


(Thus, z — ► 1 as the dibaryon density goes to zero.) This specifies the dibaryon equation 
of state in parametrized form. (Alternatively, eliminating z , the equation of state can 
be written 13 p H = (W9)[(l + 3p H /4p 0 ) 1/J - l] 3 -) We can then calculate the dibaryon 

chemical potential: 

= = (is) 

ng 

and the adiabatic index 


T b = 


Ph^ + pg 

PH 


( _ 4z 

\dpg / 3z — 1 


(16) 
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At low density ( n B -* 0). Th -» 2, and at high density (p > po ) I* -+ 4/3; the latter 
reflects the fact that the dibaryon EOS reduces to that of a free relativistic gas in the 
high density limit. Since Iff < r n , dibaryons generally soften the equation of state. 

Imposing chemical equilibrium between neutrons and dibaryons, we find that the 
threshold baryon number density for the appearance of dibaryons is: 


n B — 


/_mff — 2 Tn g_V /0 fm ~ 3 
\1198.88MeV / 


= 0.45 fin" 3 


/ m g - 2m n \ 1/a 
V 2,ttl \ 2 m n / 


(17) 


This threshold is independent of the self-coupling A. For mg - 2m A , the threshold baryon 
density for the appearance of H-particles is about three times nuclear matter density. 

For densities above nj,, the equation of state will be that of a neutron-dibaryon mix- 


ture: 


n B = n n + 2 rig, 


P = Pn + Pa, 

P =Pn +PH = P(«b)j 

PB = (P + P)/ n B = Mb(p)* (- 1 ®) 

Imposing chemical equilibrium, the system can be described with one parameter, for 
example the neutron density n n . 

We have constructed models of neutron stars by integrating the Oppenheimer-Volkoff 
equation of hydrostatic equilibrium using the above equation of state. The efFect of the 
appearance of dibaryons can be seen in Fig. 1, which shows the neutron star mass M vs. 
radius R for models with (solid curve) and without dibaryons (dashed curve). For this 
choice of parameters, m s = 2200 MeV and A = 8000, the maximum neutron star mass is 
reduced when dibaryons are included, a consequence of the softening of the equation of 
state. 

In Fig. 2, we plot the neutron star mass threshold for dibaryons to be present: a 
neutron star with BJ1 equation of state and a mass above the threshold will have an 
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admixture of dibaryons in its core. Note that dibaryons are present in stars of mass 
IAAMq for a wide range of H-particle masses. In particular, even if the H is unstable 

(mj> 2mA), it seems to be welcome. 

The dependence of the neutron star maximum mass on the dibaryon mass and coupling 
constant is shown in Fig. 3. The constraint that the maximum neutron star mass exceed 
1.44M©, the mass of the binary pulsar 14 PSR1913+16, implies a lower bound on the 
strength of dibaryon interactions (A) for a given dibaryon mass. In Fig. 4, we show this 
constraint on the dibaryon parameter space for the neutron-dibaryon model. Interestingly, 
this lower bound on A as a function of mg is close to the bag model estimate of A, Eqn.(6). 
We also note that if m H < 2 m n = 1879 MeV, then “dibaryon matter” would be absolutely 
stable at zero pressure; in this case, dibaryons would form ‘boson stars’ 13 (perhaps with 
a small crust of nuclear matter) with maximum mass M maa = 0.22(A/A ) ^ Tn^/mg — 
1.44Af©(A/2.8 x 10 3 ) 1 / 3 (2m n /m H ) 3 . 

Aside from their effects on neutron star structure, dibaryons may also have impor- 
tant consequences for the transport and cooling properties of neutron stars. Since they 
are a non-ideal Bose gas, the low-lying excitations will have a phonon spectrum, and 
the dibaryons will presumably act as a superfluid. In fact, one can demonstrate this di- 
rectly by considering the spectrum of small fluctuations 1 ' about the solution (7). Also, 
since they are not Pauli-blocked, neutrino emission from H-H and H-n scattering may be 
considerably enhanced over the modified URCA process, thereby accelerating the early 
stages of neutron star cooling. We conclude that dibaryons could play an important role 
in neutron star physics. 

Part of this work was completed while two of us (J.F. and A.O.) were visitors at 
the Copernicus Astronomical Center in May 1990; we thank the Center for its warm 
hospitality. J.F. and A.O. acknowledge support from the DOE and NASA (through grant 
NAGW-1340) at Fermilab, and A.O. from NSF (Grant AST-22595) and DOE (grant 
DE-FG02-90ER40606) at the University of Chicago. 
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FIGURE CAPTIONS 


Figure 1. Gravitational mass M versus stellar radius R for BJl equation of state 
(solid curve) and for BJl with dibaryons of mass m B = 2200AfeV and self-coupling 

A = 8000 (dashed curve). 

Figure 2. Threshold neutron star mass above which dibaryons appear, as a function 
of dibaryon mass. 

Figure 3. Maximum neutron star mass as a function of A for m H = 2000, 2100, 2200, 
2300 MeV. The horizontal line at 1.44M© denotes the mass of the binary pulsar. 

Figure 4. Minimum value of A required for stability of a 1.44M© neutron star, 
A mi „/1000, as a function of H-mass m B . 
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